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Electroweak baryogenesis is an interesting theoretical scenario, which de- 
mands physics beyond the Standard Model at energy scales of the order of 
Q^l the weak boson masses. It has been recently emphasized that, in the pres- 

p ■ ence of light stops, the electroweak phase transition can be strongly first 

1^ ■ order, opening the window for electroweak baryogenesis in the MSSM. For 

the realization of this scenario, the Higgs boson must be light, at the reach 
of the LEP2 collider. In this article, we compute the baryon asymmetry 
/\ ' assuming the presence of non-trivial CP violating phases in the parameters 

c^ ' associated with the left-right stop mixing term and the Higgsino mass /i. 

We conclude that a phase |sin0^| > 0.01 and Higgsino and gaugino mass 
parameters |/i| ~ M2, and of the order of the electroweak scale, are necessary 
in order to generate the observed baryon asymmetry. 
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1. The origin of the observed baryon asymmetry is one of the most fundamental open 
questions in particle physics. The Standard Model (SM) fulfills all the requirements [0] 
for a successful generation of baryon number 0, due to the presence of anomalous pro- 
cesses , which also put relevant constraints in models in which the baryon asymmetry 
is generated at a very high energy scale 0]. However, the electroweak phase transition 
is too weakly first order to assure the preservation of the generated baryon asymmetry 
at the electroweak phase transition |^, as perturbative |^, 0] and non-perturbative [^ 
analyses have shown. On the other hand, CP-violating processes are suppressed by 
powers oimf/M\Y, where nif are the light-quark masses. These suppression factors are 



sufficiently strong to severely restrict the possible baryon number generation ||^, 10 
Therefore, if the baryon asymmetry is generated at the electroweak phase transition, it 
will require the presence of new physics at the electroweak scale. 

Low energy supersymmetry is a well motivated possibility, and it is hence highly 
interesting to test under which conditions there is room for electroweak baryogene- 



sis in this scenario [|1T| , p!2| , |T3|. It was recently shown |TJ] that the phase transition 
can be sufficiently strongly first order only in a restricted region of parameter space: 
The lightest stop must be lighter than the top quark, the ratio of vacuum expectation 
values tan/? < 3, while the lightest Higgs must be at the reach of LEP2 Q. Similar re- 
sults were independently obtained by the authors of Ref. |T^. These results have been 



confirmed by explicit sphaleron calculations in the Minimal Supersymmetric Standard 
Model (MSSM) [0]. 

On the other hand, the Minimal Supersymmetric Standard Model contains, on top 
of the Cabbibo-Kobayashi-Maskawa (CKM) matrix phase, additional sources of CP- 
violation and can account for the observed baryon asymmetry. New CP-violating phases 
can arise from the soft supersymmetry breaking parameters associated with the stop 
mixing angle. Large values of the mixing angle are, however, strongly restricted in order 
to preserve a sufficiently strong first order electroweak phase transition |TT, |T2|, 13, H]. 



Therefore, an acceptable baryon asymmetry may only be generated through a delicate 
balance between the values of the different soft supersymmetry breaking parameters 
contributing to the stop mixing parameter, and their associated CP-violating phases. 
The value of the Higgsino and gaugino masses also play a very relevant role in the 
determination of the baryon asymmetry. Indeed, the final baryon asymmetry depends 
strongly on the relative value of the Higgsino and gaugino mass parameters with respect 
to the critical temperature Tc = C(100 GeV). 

In this letter we compute the baryon asymmetry and the strength of the first order 
phase transition in the MSSM. We identify the region in the supersymmetric parameter 
space where baryon asymmetry is consistent with the observed value ns/s ~ 10^^'^ 
and, furthermore, it is not washed out inside the bubbles after the phase transition. We 
obtain that the second effect is guaranteed provided the light stop mass is in the range 
Mz < mj ^ rrit, the lightest Higgs boson mass is bounded by mn ~ 80 GeV and the 
CP-odd boson has a mass m^ > 150 GeV, independently of the chargino and neutralino 
masses. On the other hand baryon asymmetry is generated mainly by charginos and 
neutralinos which are not much heavier than the critical temperature, independently of 



^Explicit two- loop calculations have the general tendency to strengthen the phase transition [|5| |T^ 
thus making the previous bounds as very conservative ones. 



the mass of the hghtest stop. 



2. As discussed above, a strongly first order electroweak phase transition can be achieved 
in the presence of a top squark hghter than the top quark |]14|. In order to naturally 



suppress its contribution to the parameter Ap and hence preserve a good agreement 
with the precision measurements at LEP, it should be mainly right handed. This can 
be achieved if the left handed stop soft supersymmetry breaking mass mq is much larger 
than Mz- For moderate mixing, the lightest stop mass is then approximately given by 



Ai 



m\c^mij + Di, + mi{<P) | 1 - ^2 



where At = At — yU*/tan/5 is the particular combination appearing in the off-diagonal 
terms of the left-right stop squared mass matrix and m\j is the soft supersymmetry 
breaking squared mass parameter of the right handed stop. Observe that the stop mass 
eigenvalues and hence, the stop contribution to the Higgs effective potential must be 
computed by taking into account the non- vanishing complex phases for the parameters 
At and /x. For most practical purposes, however, unless the phases are of order one, we 
can identify the absolute value of the parameters At and /i with their real components 
and hence, we shall not distinguish them in our notation. 

We shall work within the framework of the improved one-loop effective potential. 
Hence two loop corrections |TB|, |TB| will not be included in our analysis. As stated 



above, since these corrections tend to increase the strength of the first order phase 
transition |jT6|, our bounds should be taken as conservative ones. The preservation 
of the baryon number asymmetry requires the order parameter v{Tc)/Tc to be larger 
than one, where v denotes the Higgs vacuum expectation value (VEV) and we have 
normalized it such that v{0) = v = 246.22 GeV. Following the analysis of Refs. [p!2| , [14 



the order parameter v{Tc)/Tc is bounded to be below the maximum value obtained for 
iriA > Tc, 




(2) 
IT V mj 

where rrit = mt{mt) is the on-shell running top quark mass in the MS scheme. The first 
term on the right hand side of Eq. (^ is the Standard Model contribution 

-^ -i '° V (3) 

TJsM \mH\GeY]) ' ^ ' 

and the second term is the contribution that would be obtained if the right handed 
stop plasma mass vanished at the critical temperature (see Eq. @)). The lower the 
value of this mass at the critical temperature, the larger the value of v{Tc)/Tc |T^. 
For rriji > 150 GeV, preferred to enhance the strength of the phase transition [jl2|, the 
present experimental bound on the Higgs mass reads mn > 65 GeV, meaning that the 
SM contribution leads to a value of v{T^/T(. at most of order 1/3. 



In order to overcome the Standard Model constraints, the stop contribution must be 
large. The stop contribution strongly depends on the value of m^, which must be small 
in magnitude, and negative, in order to induce a sufficiently strong ffist order phase 
transition. Indeed, large stop contributions are always associated with small values of 
the right handed stop plasma mass 

mf = -ml + I[niT) (4) 

where m^ = -mfj, IIr{T) ~ Ag^T^ /9 + h^ 6[2 - A^Jml]T^ MS is the finite temper- 
ature self-energy contribution to the right-handed squarks and ht and gs are the top 
quark Yukawa and strong gauge couplings, respectively. Moreover, the trilinear mass 
term. At, must be A^ -C rriQ in order to avoid the suppression of the stop contribution 
to v{Tc)/Tc. 

Although large values of niu, of order of the critical temperature, are useful to get 
a strongly ffist order phase transition, they may also induce charge and color breaking 
minima. Indeed, if the effective plasma mass at the critical temperature vanished, the 
universe would be driven to a charge and color breaking minimum at T > T^ |T| 



Hence, the upper bound on v{Tc)/Tc, Eq. (||) cannot be reached in realistic scenarios. 
A conservative bound on niij may be obtained by demanding that the electroweak 
symmetry breaking minimum should be lower than any color-breaking minima induced 
by the presence of rriu at zero temperature, which yields the condition 

rnu<{^y'. (5) 

It can be shown that this condition is sufficient to prevent dangerous color breaking 



minima at zero and finite temperature for any value of the mixing parameter At |T4 
In this work, we shall use this conservative bound. 

In order to obtain values of v{Tc)/Tc larger than one, the Higgs mass must take 
small values, close to the present experimental bound. Numerically, an upper bound, 
of order 80 GeV, can be derived. For small mixing, the one-loop Higgs mass has a very 
simple form 



ml = M|cos2 2/3 + -^^log -^^ 



l + O 




(6) 



where m~ ~ niQ + mf, is the heaviest stop squared mass. Hence, tan/3 must take values 
close to one. The larger the left handed stop mass, the closer to one tan/5 must be. In 
this work we shall take mq > 500 GeV. This implies that the left handed stop effects 
decouple at the critical temperature and hence, different values of tuq mainly affect the 
baryon asymmetry through the resulting Higgs mass. 



3. Concerning the baryon asymmetry generation, the new source of CP- violation, be- 
yond the one contained in the CKM matrix, may be either explicit ||2^ or sponta- 
neous l^ll in the Higgs sector (which requires at least two Higgs doublets). In both 



^We are considering heavy (decoupled from the thermal bath) gluinos. For light gluinos, their 
contribution to the squark self-energies, 2g^T^/9, should be added to Ilii{T) 113]. 



cases, particle mass matrices acquire a nontrivial space-time dependence when bubbles 
of the broken phase nucleate and expand during a first-order electroweak phase tran- 
sition. The crucial observation is that this space-time dependence cannot be rotated 
away at two adjacent points by the same unitary transformation. This provides suf- 
ficiently fast nonequilibrium CP-violating effects inside the wall of a bubble of broken 
phase expanding in the plasma and may give rise to a nonvanishing baryon asymmetry 
through the anomalous {B + L)-violating transitions [Q when particles diffuse to the 
exterior of the advancing bubble. 

Baryogenesis is therefore fueled by CP-violating sources which are locally induced by 



the passage of the wall [23, "M]. These sources should be inserted into a set of classical 
Boltzmann equations describing particle distribution densities and permitting to take 
into account Debye screening of induced gauge charges |2^ 0, particle number changing 
reactions |2^ and to trace the crucial role played by diffusion [^. Indeed, transport 
effects allow CP-violating charges to efficiently diffuse in front of the advancing bubble 
wall where anomalous electroweak baryon violating processes are unsuppressed. This 
amounts to greatly enhancing the final baryon asymmetry. 

We shall make use of a method recently proposed by one of the authors [Q to com- 
pute the effect of CP-violation coming from extensions of the Standard Model on the 
mechanism of electroweak baryogenesis. This method is entirely based on a nonequilib- 
rium quantum field theory diagrammatic approach and has two main virtues: it may 
be applied for all wall shapes and sizes of the bubble wall and it naturally incorporates 
the effects of the incoherent nature of plasma physics on CP-violating observables. The 
method relies on the closed-time path (CPT) formalism, which is a powerful Green's 



function formulation for describing nonequilibrium phenomena in field theory ||2^, [30|] , 
as the ones we are interested to occur in the neighbourhood of the advancing bubble 
walls. Indeed, what we need is to compute the temporal evolution of classical order 
parameters, namely CP-violating particle currents, with definite initial conditions. In 
this respect, the ordinary equilibrium quantum field theory at finite temperature may 
not be applied, since it mainly deals with transition amplitudes in particle reactions. 

It is well-known that self-energy corrections at one- or two-loops to the propagator 
modify the dispersion relations and introduce nontrivial effects {e.g. damping) due to 
the imaginary contributions to the self-energy pT| . This is because particles propagating 



in the plasma experience interactions with the surrounding particles of the thermal 
bath, and the correct way to describe this phenomenon is to substitute particles by 
quasiparticles and to adopt dressed propagators. The self-energy takes the form S(fc) = 
Re S(/c) -|- i Im S(A;). Due to the nonvanishing Im S, the spectral function p(k, /c°) 
acquires in the weak limit a finite width 

2 uj[k) 
uj\k) = k2 + m2 + ReS(k,cu), (7) 



^This effect will be ignored since the impact on the final result is 0(1) [^ 



and is expressed by 
p(k, k') = 



{ko + ie + iVy- uj^ik) {k^-ie-iTy- iu^{k) 



where m indicates the tree-level mass. It is easy to show that the free propagator is 
obtained when taking the limit F ^ 0. 

The advantage of the method is that very general formulae can be obtained for 
the temporal evolution of CP- violating observables without any particular assumption 
on the relative magnitude of the mean free paths and the thickness of the wall. The 
standard thin wall and thick wall limits are recovered in the limit TL^^ — > and oo, 
respectively. Moreover, because of the presence of a finite damping rate, the Green's 
functions of the particles involved in the processes are damped for times t > T~^, 
reflecting the effects of the incoherent nature of the plasma on CP- violating observables. 

Following [^, |3^, we are interested in the generation of charges which are approx- 
imately conserved in the symmetric phase, so that they can efficiently diffuse in front 
of the bubble where baryon number violation is fast, and non-orthogonal to baryon 
number, so that the generation of a non-zero baryon charge is energetically favoured. 
Charges with these characteristics are the axial stop charge and the Higgsino charge, 
which may be produced from the interactions of squarks and charginos and/or neutrali- 
nos with the bubble wall, provided a source of CP-violation is present in these sectors. 
CP- violating sources 7q(-z) (per unit volume and unit time) of a generic charge density 
J° associated with the current J'^(z) and accumulated by the moving wall at a point 
z'^ of the plasma can then be constructed from J'^^z) |^ ^ 

7q(z) = doJ%z). (9) 

The definition of the CP- violating source •yqiz) is appropriate to describe the damping 
effects originated by the plasma interactions, but does not incorporate any relaxation 
time scale arising when diffusion and particle changing interactions are included. How- 
ever, one can leave aside diffusion and particle changing interactions and account for 
them independently in the rate equations ||32|, |33| . 

3.1 The right-handed stop current J^ associated to the right-handed stop t/j is given 

r^ = z {t},d^tn - d^t*^tn) . (10) 

{Jr(z)) gets contributions from eight different one-loop triangle Feynman diagrams, 
symbolically depicted in Fig. la. 

In the calculation of {J{i{z)) at the lowest order in a "Higgs insertion expansion", 
and more precisely in the integrand of / dxdy {x and y being the time components of 
the four- vectors where the two insertions of the diagram in Fig. la have been defined 0), 
the function H2{x + z)Hi{y + z) — Hi{x + z)H2{y + z) appears, where Hi, i = 1,2, 
denote the two neutral Higgs components. In order to deal with analytic expressions, 

^Hereafter the generic notation for the time component of the four- vector z^ will be adopted as z, 
e.g. z^ = (z,z). 



we can work out the thick wall limit and simplify the expressions obtained above by 
performing a derivative expansion 

oo 1 Qn 

^^(^ + -) = E-r53r^^(^)^"- (11) 

n=0 "'■ ^^ 

The contribution to the current can then be expanded in a power series as: 

H2ix + z)Hi{y + z)- Hi{x + z)H2{jj + z) 

oo On 

= Y.^^[Hi{z)9zH2{z)-H2{z)d.H,{z)]Ux,y) (12) 

n=0 '^^ 

where the function /„ is provided by the power expansion, e.g. fo = x — y, fi = 
(x^ — y^)/2, . . . Notice that the term with no derivatives vanishes in the expansion (|T^, 
e.g. H2{z)Hi{z) — Hi{z)H2{z) = 0, which means that the static term in the derivative 
expansion (pUf ) does not contribute to the vacuum expectation value of the currents. 



Furthermore, the approximation of neglecting terms in (O) with n > 1 amounts to 



neglecting terms in ([T2|) with n > 0. Since for a smooth Higgs profile, the subsequent 
derivatives with respect to the time coordinate are associated with higher powers of 
Vw/Luj, and the integration over x and y of a higher order term in the expansion of the 
currents leads to a higher power of l/F^, this is a good approximation for values of 
L^Tj/vi^ ^1. In other words, this expansion is valid only when the mean free path 
Tj ~ r~^ is smaller than the scale of variation of the Higgs background determined by 
the wall thickness, L^, and the wall velocity f^j. An estimate of the stop damping rate 
(in the low momentum limit) is r^~ 10"^ T, approximated using the calculation made 



in Ref. [^. With such value, our derivative expansion is perfectly justified since the 
wall thickness can span the range (10 — 100)/T. 

Therefore, the currents are proportional to the function (the coefficient of /o in Eq. 

(0)) 

H,{z)d,H2{z) - H2{z)d,H,{z) = H\z)d,(3{z), (13) 

which should vanish smoothly for values of z outside the bubble wall. Here H^ = 
Hi + Hi 

Since the time variation of the Higgs fields is due to the expansion of the bubble 
wall through the thermal bath, {Jji{z)) will be linear in f^. This result explicitly shows 
that we need out of equilibrium conditions to generate {J^iz)) and that we have to 
call for the CTP formalism to deal with time- dependent phenomena. To work out 
exactly {J'^{z)) one should know the exact form of the distribution functions which, 
in ultimate analysis, are provided by solving the Boltzmann equations. However, any 
departure from thermal equilibrium distribution functions is caused at a given point by 
the passage of the wall and, therefore, is 0{v^). Since {Jr{z)) is already linear in v^, 
working with thermal equilibrium distribution functions amounts to ignoring terms of 
higher order in f^ |Q, which is as accurate as the bubble wall is moving slowly in the 
plasma and we shall adopt this approximation from now on. 

One can show that 

(J°(^)) = h^,lm{Atfi)[Hi{z)d,H2{z)-H2{z)d,Hi{z)]g'^ (14) 



where 



Q 



Gl 



4:7l'^Tii Uq Ur 



(1 + 2Re(nQ)) /i(cu^, Tn, ujq, Tq) + (1 - 
- 2 {lm{nji) + lm{nQ)) hiujR, Tr, uq, Tq) 



2Re{nR)) h{ujQ, Tq, ujr, Tr) 



and a;|(^) 



^^ + ^u{Q) + nc/(Q), n^(Q) = 1/ exp (^^(Q)/r + ir^,(Q)/T) - 1 



i? ~ J- Q ~ J- 7- 
The functions Ji and I2 are given by 



/i(a, 6, c, d) 



'"I 



+ 



r2 



{rl + 1) [(a + c)2 + (6 + d)2] (ri + 1) [(a - cf + {b + d) 



21' 



(15) 



while 



(16) 



l2{a,b,c,d) 



+ 



ri-1 



2 (rf + 1) [(a + c)2 + (6 + d)2] 2 (ri + 1) [{a - cY + (6 + dy] 
where ri = (a + c)/(6 + d) and r2 = (a — c)/(6 + d). 



(17) 



3.2 Similarly, the Higgsino current associated with neutral and charged Higgsinos 

can be written as 

~ (18) 






where H is the Dirac spinor 



H 



£2 
Hi 



(19) 



and H2 = H2 {^2)1 Hi = H^ {Hi) for neutral (charged) Higgsinos. 

The Higgsino current ([I8D gets its contribution from the typical triangle diagram of 
Fig. lb and may be computed along the same lines described for the axial stop number 
in the previous subsection. Analogously to the case of right-handed stops, the dispersion 
relations of charginos and neutralinos are changed by high temperature corrections ||35|| . 
Even though fermionic dispersion relations are highly nontrivial, relatively simple ex- 
pressions for the fermionic spectral functions may be given in the limit in which the 
damping rate is smaller than the typical self-energy of the fermionic excitation (this limit 
is certainly satisfied in our case) [^ . For instance, the spectral function of Higgsinos 
H (and analogously for gauginos W, B) may be written as 



P5(k,fc° 



(A + 



m 



H 



(20) 



[k^ + ie + i 1^)2 - cu|(A;) {k^-ie-i 1^)2 - u\{k) 



where oj'~^{k) = k^ + m~(T) and m?~{T) is the Higgsino effective plasma squared mass in 
the thermal bath which may be well approximated by its value in the present vacuum. 



m'~{T) ^ |yup. Similarly, \fi\ should be replaced by M2 for p^{k,k^), and by Mi for 
P5(k,A;0). 

In analogy to the stop case, in the region of parameters defined above, we can perform 
a Higgs insertion expansion of the CP- violating current. The vacuum expectation value 
of the (zero component of the) Higgsino current is then found to be 



{J'~{z)) = Im(^) [H,{z)d,H2{z) - H2{z)d,Hi{z) ] ^3M2 gl ^f + Mi g^ gS 
where 



(21) 



-,w 



00 j^2 

dk — 

Ati^'TjjU^u^ 



1 - 2Re(n^)) h{ujj^, T^, w^, r^;^;.) + (l - 2Re(n^)) h{uj^, T 



W^H'> H> 



+ 2 (im(n^) + Im(n^)) hi^uj^, V 



H'^W W' 



^^^ ^H(w) = ^^ + l/^P(^l) while nj^^^^^ = 1/ [exp (u^^^^^/T + iV^^^^lT) + 1 



:22i 



. The 



H(W) ■- ' ""I ^^■^''■' --HiW) 

damping rate of charged and neutral Higgsinos is dominated by weak interaction and 
we take F^ ~ F^ to be of order of 5 x 10~^T. The Bino contribution may be obtained 
from the above expressions by replacing M2 by Mi. 



4. We can now solve the set of coupled differential equations describing the effects 
of diffusion, particle number changing reactions and CP-violating source terms. We 
will closely follow the approach taken in Ref. |23| where the reader is referred to for 
more details. If the system is near thermal equilibrium and particles interact weakly, 
the particle number densities Ui may be expressed as rii = kifiiT^/6 where /Xj is the 
local chemical potential, and ki are statistical factors of order of 2 (1) for light bosons 
(fermions) in thermal equilibrium, and Boltzmann suppressed for particles heavier than 
T. 

The particle densities we need to include are the left-handed top doublet ql = 
[ih + &l), the right-handed top quark tn, the Higgs particle h = (if°, H2, -fff , -^2*^), and 
the superpartners qi, Ir and h. The interactions able to change the particle numbers 
are the top Yukawa interaction with rate F^, the top quark mass interaction with rate 
Fm, the Higgs self-interactions in the broken phase with rate F-^, the strong sphaleron 
interactions with rate Fgs, the weak anomalous interactions with rate F^s and the gauge 
interactions. We shall assume that the supergauge interactions are in equilibrium. 
Under these assumptions the system may be described by the densities Q = Ql + Ql, 
T = tn + tji and 7i = h + h. CP-violating interactions with the advancing bubble 
wall produce source terms 7^ = dz{J~{z)) for Higgsinos and jr = dz{JR{z)) for right- 
handed stops, which tend to push the system out of equilibrium. Ignoring the curvature 
of the bubble wall, any quantity becomes a function of the coordinate z = Z3 + v^^z, 
the coordinate normal to the wall surface, where we assume the bubble wall is moving 
along the z^-axis. 



8 



Assuming that the rates Tf and Fss are fast so that Q/kq — Ti/k-^ — T/kq- = 0{l/rt) 
and 2Q/kg — T/kr + 9(Q + T)/kb = 0{1/Tss), one can find the equation governing the 
Higgs density 

v^n' - DH" + TTi: - 7 = 0, (23) 

where the derivatives are now with respect to z, D is the effective diffusion constant, 
7 is an effective source term in the frame of the bubble wall and F is the effective 
decay constant [^. An analytical solution to Eq. (|23|) satisfying the boundary con- 



ditions 7i(±oo) = may be found in the symmetric phase (defined by z < 0) using 
a z-independent effective diffusion constant (in agreement with our approximation in 
Eq. ([TT|) ) and a step function for the effective decay rate F = r^(z). A more realistic 
form of F would interpolate smoothly between the symmetric and the broken phase 
values. We have checked, however, that the result is insensitive to the specific position 
of the step function inside the bubble wall. The values of D and F in (^) of course de- 
pend on the particular values of supersymmetric parameters. For the considered range 
we typically find 13 ~ 0.8 GeV'S V ~ 1.7 GeV. 

The dependence on 7 is much more subtle, and a good approximation for the 
Higgs profiles is necessary in order to obtain a reliable approximation to the cur- 



rents. As seen in Eq. (|T3|), the current is proportional to the variation of the ratio 
of vacuum expectation values along the wall. The Higgs profiles along the wall are 
likely to follow the path of minimal energy connecting the false and the true vacua 
in the Higgs potential. Close to the symmetric phase, and at the temperature Tq, 
defined as the one at which the curvature at the origin vanishes, the ratio of vac- 
uum expectation values is given by tan^/5(z = 0) = m\{TQ) / m\{TQ) , a quantity which 
may be computed numerically from the Higgs finite temperature effective potential as 
m?i{T) = 1/2 (d'^Vcs/dHf) 1^^=^2=05 (^ = 1,2). Hence, the variation of (3 is given by 

AP = /?(To) - arctan(mi(To)/m2(ro)) (24) 

This quantity tends to zero for large values of rriA^, and takes small values, of order 0.015 
for values of rriA = 150-200 GeV, giving rise to a further suppression of the final baryon 
number. 

One would be tempted to adopt a linear ansatz already for the function (|13]) in the 
current, as done in [^. However, this approximation has the drawback that it may 
not be obtained starting from smooth Higgs profiles and can therefore be viewed as 
unphysical. Also, since the source in the diffusion equation depends on the derivatives 
of the current, linear walls (or, generally, all non-differentiable wall shapes) lead to 
the appearence of (5-functions in the CP-violating source. A more realistic option is to 
specify the shapes of the Higgs fields and /3 and to consequently compute the source. 
Realistic shapes for H{z) (and probably for /3(z)) would be given by functional kinks 
along the bubble wall. This provides a smooth function for the source 7. Moreover, 
this way of proceeding leads to a strong suppression with respect to the result that one 
would have obtained assuming a linear approximation for the quantity Hi{z)dzH2{z) — 
H2{z)dzHi{z). This is due to the fact that, since Hi{z)dzH2{z) — H2{z)dzHi{z) goes to 
zero for z — *■ ±00, the source 7 must change sign inside the bubble wall, giving rise to 
some cancellation for the final expression of the baryon number, see Eqs. (29) and (32). 



We consider here the following semirealistic approximation: 



(3{z) 



1 — cos 



zvr 



[e{z) 



Af3 



cos 



zvr 



[9iz) 



-e{z-L^)] + v{L^)e{z-L^) 

z - L^)] + f3iz = 0) + A/3e{z - L^) (25) 



which leads to a current vanishing smoothly at z = and z 
For z < one obtains, 



n{z) =^e^^-/^, 



(26) 



and for z > 0, 



n{z) 



B^-^ 



1 



+ S- 



D{X+ - A_) ^0 
D{\_ - A+) io 



'du7(u)e~^-" 



^A+z 



.A_z 



where 



A. 



2D 



(27) 



(2^ 



and 7(z) = Vi^dzJ'^{z)f{ki), Jq being the total CP-violating current resulting from 
the sum of the right-handed stop and Higgsino contributions and /(fcj) a coefficient 
depending on the number of degrees of freedom present in the thermal bath and related 



to the definition of the effective source |3^. Imposing the continuity of H and H' at 
the boundaries, we find 



A = B, 



A+ 



i3_ 



'A+ 



1 



1 =^ 



D A_ 



du ■y{u)e 



'\+U 



(29) 



From the form of the above equations one can see that CP-violating densities are non 
zero for a time t ~ D/v'^ and the assumptions leading to the analytical form of 7i(z) 
are valid provided Fj, Fgg ^ v'^/D. 

The equation governing the baryon asymmetry ub is given by ||33| 



DqfiB - v^n'j^ 



-z)n/Fws^L = 0, 



(30) 



where F^s = 6fi;a^T is the weak sphaleron rate (k ~ 1) [^305 ^^^ rii is the total number 
density of left-handed weak doublet fermions, rij = 3 is the number of families and we 
have assumed that the baryon asymmetry gets produced only in the symmetric phase. 
Expressing ^^(z) in terms of the Higgs number density 



ul 



9kqkT — Skbkq — 5kbkq 
knih + Qkq + 9kr) 



n 



(31) 



^We have checked that numerical results using the ansatz ( p5| ) differ by less than 10% from those 
obtained by using functional kinks. This leads us to the belief that our results are robust even in the 
absence of a numerical calculation of bubble solutions in the MSSM. 

^The correct value of k is at present the subject of debate; see, for instance, Ref. |3q|. 
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and making use of Eqs. (p6D-(pOD, we find that 

— = -9{ki) — 5 , (32) 

s vis 

wliere s = 27i'^g^,sT^/4:5 is the entropy density {g^s being the effective number of rela- 
tivistic degrees of freedom) and g{ki) is a numerical coefficient depending upon the hght 
degrees of freedom present in the thermal bath. 



5. The expression for the baryon asymmetry derived in the last section assumes that 
the order parameter v{Tc)/Tc > 1, in order to assure that anomalous processes are 
sufficiently suppressed in the broken phase. As we already discussed, this demands low 
values of the right-handed stop effective plasma mass at the critical temperature and 
low values of the mixing mass parameter At. We will now present numerical results 
based on the previous calculations. Since A/5, Eq. (P^, denotes the variation of the 
ratio of vacuum expectation values of the Higgs fields along the bubble wall, it will 
be negligible unless the CP-odd Higgs mass takes values of the order of the critical 
temperature. Indeed, for large values of niA, A/? ~ ^/fn\. 

Values of the CP-odd Higgs mass rriA ~ 200 GeV are associated with a weaker first 
order phase transition. Fig. 2 shows the behaviour of the order parameter v/T in the 
m^-tan/3 plane, for At = 0, ttiq = 500 GeV and values of mj/ close to its upper bound, 
Eq. (^. In order to correctly interpret the results of Fig. 2 one should remember that 
the Higgs mass bounds are somewhat weaker for values of rriA < 150 GeV. However, 
even for values of ttia of order 80 GeV, in the low tan /3 regime the lower bound on the 
Higgs mass is of order 60 GeV. Hence, it follows from Fig. 2 that, to obtain a sufficiently 
strong first order phase transition, v/T > 1, the CP-odd Higgs mass tua ~ 150 GeV f\. 
In our analysis, we shall fix rriA = 200 GeV, which is well inside the acceptable range 
for this parameter. As anticipated above, this choice leads to values of A/3 ^ 0.015. 

Since for large values of rriQ the stop contribution to the baryon asymmetry is 
strongly suppressed compared to the chargino and neutralino ones, the numerical values 
of the baryon asymmetry depend linearly on the phase of the Higgsino mass parameter. 
In general, the numerical values obtained following the above procedure are too low 
to explain the observed baryon asymmetry unless this phase is of order one. This 
conclusion may only be avoided for certain specific values of the gaugino and Higgsino 
mass parameters. In order to make this result explicit, in Fig. 3 we show the value of this 
phase necessary to obtain a value ns/s ~ 4 x 10~^^ in the H-M2 plane. The wall velocity 
is taken to be f^^ = 0.1, while the bubble wall width is taken to be L^^ = 25/T, and, for 
simplicity, we consider Mi = M2 [|. Our results, are, however, quite insensitive to the 
specific choice of v^ and L^. Due to the fact that, as already stated, stop contributions 
to the baryon asymmetry are negligible as compared to the chargino/neutralino ones, 
we have taken | sin(<^^ + 0„)| = in Fig. 3. 



"^Recent analyses [39| of the effective three-dimensional theory of the MSSM essentially confirm the 
early results of Ref. |fj], while a similar analysis for the light stop scenario is still lacking. 

^We have checked that, as expected, the final results do not depend much on the value of Mi. In 
particular Mi can be chosen such that the lightest neutralino is lighter than the light stop. 
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Values of the phases lower than 0.1 are only consistent with the observed baryon 
asymmetry for values of |/i| of order of the gaugino mass parameters. This is due 
to a large enhancement of the computed baryon asymmetry for these values of the 
parameters. Since the Higgs background in the Feynman diagrams of Fig. lb is carrying 
a very low momentum (of order of the inverse of the bubble wall width L^^), this resonant 
behaviour is associated with the possibility of absortion (or emission) of Higgs quanta 
by the propagating particles. For momenta of order of the critical temperature, this 
can only take place when the Higgsinos and gauginos are nearly degenerate in mass, 
/i ~ M2. By using the Uncertainty Principle, it is easy to understand that the width of 
this resonance is expected to be proportional to the decay width of the -inos, which is 
what comes out from our explicit computation. In the case of the stops, instead, since 
ttiq ^ Tnfj, the same phenomenon can only happen for momenta larger than mq. Such 
configurations are exponentially suppressed and do not give any relevant contribution 
to the one loop diagram of Fig. la. The resonant behaviour is certainly enhanced by our 
approximations. Indeed, since we are working at first order in the Higgs insertions, the 
relevant masses propagating in the loops are given by their expressions in the symmetric 
phase. Hence, strong degeneracies are possible when /i is close to M2. This degeneracy, 
however, will be broken by further Higgs insertions. 

Therefore, taking the results from Fig. 3 at face value, one would find an absolute 
lower bound on |sin(0^)| > 10~^. However, given all uncertainties in the calculation, 
including those associated to the different damping rates for squarks, Higgsinos and 
gauginos, and the sphaleron parameters, in particular k = T^s/[6af^T] ^ we will prefer 
to consider, as a conservative estimate, the previous bound on | sin(0^)| accurate only 
up to about one order of magnitude. If the phases of the parameters At and n are 
0{1), CP- violation will be induced in the Higgs sector and the Higgs VEVs will become 
complex. However, it is easy to prove that, whenever the tree level CP-odd mass is 
larger than the vacuum expectation values and the trilinear mass parameters l^^l and 
|/i| take values lower than mq, the induced phase of the product of vacuum expectation 
values will be much lower than 1, hence our results will be only weakly modified by this 
effect. 



6. In conclusion, we made a self-consistent computation of the baryon asymmetry gen- 
erated at the electroweak phase transition in the Minimal Supersymmetric Standard 
Model. We proved that baryon asymmetry is mainly generated by CP-violating Hig- 
gsino currents, provided that Higgsinos and gauginos are not much heavier than the 
electroweak critical temperature {Tc ~ 100 GeV), and the phase 0^ is not much smaller 
than 0.1. Observe that these relatively large values of the phases are only consistent 
with the constraints from the electric dipole moment of the neutron if the squarks of 
the first and second generation have masses of the order of a few TeV |^. On the 



other hand, the generated baryon asymmetry is not washed out inside the bubbles, 
after the phase transition, provided that the lightest stop is lighter than the top quark 
(the so-called light stop scenario), the pseudoscalar Higgs boson is heavier than ~ 150 
GeV and the lightest Higgs boson is lighter than ~ 80 GeV. 

® We have adopted in the calculation the value k = 1. The final result scales with k. 
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The most direct experimental way of testing this scenario is through the search 
for the hghtest Higgs at LEP2. If the Higgs is found, the second test will come from 
the search for the lightest stop at the Tevatron collider (the stop mass is typically 
too large for the stop to be seen at LEP). If both particles are found, the last crucial 
test will come from B physics, more specifically, in relation to the CP-violating effects. 
Depending on their gaugino-Higgsino composition, light neutralinos or charginos, at 
the reach of LEP2, although not essential, would further restrict the allowed parameter 
space consistent with electroweak baryogenesis. 

Finally, let us expand on the Higgs searches. If the next run of LEP at ^/s = 186 
GeV collects the planned luminosity, it will be able to search for Higgs masses up to 
about 75 GeV [0]. If no Higgs signal is found, this will pose additional constraints on 
the present scenario. In particular, in order to get still a sufficiently strong first order 
phase transition it will demand At/rriQ ~ 0. Moreover, the selected parameter space 
leads to values of the branching ratio BR(6 — > sj) different from the Standard Model 
case. Although the exact value of this branching ratio depends strongly on the value 
of the /i and At parameters, the typical difference with respect to the Standard Model 
prediction is of the order of the present experimental sensitivity and hence in principle 
testable in the near future. Indeed, for the typical spectrum considered here, due to the 
light charged Higgs, the branching ratio BR(6 — >■ 57) is somewhat higher than in the 
SM case, unless negative values of A^fi are present. 
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Figure 1: a) One- loop Feynman diagrams contributing to {Jfiiz)). Here Tn\^ indicates 
the combination ht {AtH2 — ^*Hi). b) One-loop diagrams contributing to {J~^{z)). Here 

, 4 describes the interactions of Higgsinos with gauginos Wa, a = 1, ■ ■ ■ , 3 

Kit 75), 



and B, a 
9a = 92, a 



4. More precisely, Ha = 9a 
1, ■ ■ ■ , 3, and ga = 91 for a = 



4. 
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Figure 2: Contour plots of constant values of v{Tc)/Tc (solid lines) and rriH in GeV 
(dashed lines) in the plane (m^,tan/9). We have fixed rrit = 175 GeV and the values of 
sypersymmetric parameters: rriQ = 500 GeV, rriu = m^'* fixed by the charge and color 
breaking constraint, and At = yu/tan/?. 
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Figure 3: Contour plot of | sin0^| in the plane (/i, M2) for fixed Ub/ s = 4 x 10 ^^ and 
Vo; = 0.1, L^^ = 25/T, rriQ = 500 GeV, mu = m™^*, tan/3 = 2 and At = /i/tan/?. 
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